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Abstract 

In this paper, we find a new large scale instability displayed by a ro¬ 
tating flow in forced turbulence. The turbulence is generated by a small 
scale external force at low Reynolds number. The theory is built on the 
rigorous asymptotic method of multi-scale development. The nonlinear 
equations for the instability are obtained at the third order of the pertur¬ 
bation theory. In this article, we explain a detailed study of the nonlinear 
stage of the instability and generation vortex kinks. 
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1 Introduction 

It is well known, that the rotating effects play an important role in many 
practical and theoretical applications for fluid mechanics [T] and are espe¬ 
cially important for geophysics and astrophysics [2], m, when one have to 
deal with rotating objects such as Earth, Jupiter, Sun etc. Rotating fluids 
could generate different wave and vortex motions. For example, gyroscopic 
waves, Rossby waves, internal waves, located vortices and coherent vortex 
structures m-m- Among the vortex structures the most interesting are 
the large scale ones, since they carry out the efficient transfer of energy 
and impulse. The structures which have characteristic scale much more 
than the scale of turbulence or of the external force which generates this 
turbulence are understood as large scale ones. At present we can note, 
that there are a lot of instabilities which generate the large scale vortex 
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structures (see for example m-M), including in rotating fluid with the 
non-homogeneous turbulence dn. In this work we find new large scale in¬ 
stability in rotating fluid, under the impact of small external force which 
keeps up turbulent fluctuations. The nonlinear big scale helical vortex 
structures of Beltrami type or localized kinks with internal helical struc¬ 
ture appear as a result of the development of this instability in rotating 
fluid. We can consider that external small scale force substitutes the ac¬ 
tion of small scale turbulence. It is supposed that external force is in 
plane {X,Y), which is perpendicular to the rotation axis, i.e. axis Z is 
directed along the vector of angular velocity of rotation fl. Helical 2D 
field of velocity W^, Wy turns around axis Z, when Z changes in the kink 
which links the hyperbolic point and the stable focus (figl2}. Moreover 
this field does some turns in the kink, which links instable and stable 
focuses (Fig|3]). The found instability belongs to the class of instabilities 
called hydrodynamic a-effect. For these instabilities the positive feedback 
between velocity components of following type is typical: 

drW^ - AW, - Oy^Wy = 0, 

drWy - XWy + = 0 , 

what leads to the instability, a-effect origins from magnetic hydrodynam¬ 
ics, where it engenders the increase of large scale magnetic fields (see for 
example [IS]). Afterwards it was generalized on usual hydrodynamics. 
Some examples of hydrodynamics a-effect [8|-|14| known for the time be¬ 
ing. From this point of view, a new example of a-effect is found in this 
work. The theory of this instability is developed rigourously using the 
method of asymptotic multi-scale development similar to what was done 
by Frisch, She and Sulem for the theory of the AKA effect [13]. This 
method allows to find the equations for large scale perturbations as the 
secular equations of perturbation theory, to calculate the Reynolds stress 
tensor and to find the instability. The small parameter of asymptotical 
development is the number of Reynolds R, R 1. Our paper is organised 
as follows: in Section [2] we formulate the problem and the main equations 
in rotating system of coordinates; in Section Owe examine the principal 
scheme of the multi-scale development and we give the secular equations. 
In Section [4] we calculate the velocity field of zero approximation. In Sec- 
tionOwe describe the calculations of the Reynolds stress and find the large 
scale instability. In Section [6] we discuss the saturation of the instability 
and find non linear stationar vortex structures. The results obtained are 
discussed in the conclusions given in Section [7] 
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2 The Main Equations and Formulation 
of the Problem 

Let us examine the equations of motion for non-compressible rotating fluid 
with external force Fq in rotating coordinates system : 

(vv)v + 2nxV ^-—VP + vAV + Fo, (1) 

at \ / po 


divV = 0 (2) 

The external force Fq is divergence-free. Here H-angular velocity of fluid 
rotation, viscosity, po-constant fluid density. Let us design characteristic 
amplitude of force /o, and its characteristic space and time scale Ao and 
to respectively. 

Then To = foFo • We will design the characteristic amplitude 

of velocity, generated by external force as vo- We choose the dimensionless 
variables (t,x,V): 


X 


t 

to ’ 


V 


Vq ’ 


Fo 


To 

fo’ 


PoPo 


, _ K r> _'^ov VoV 

^0 ? -^0 ^ ? Jo \ 2 ’ 

V Ao Ag 


Vo 


fo\l 


Then, in dimensionless variables the equation m takes forme : 


(j V/ —# —# —» —» _ —» —» 

^ + R{V-V)V + DxV = -VP + AV + Fo, (3) 

g _ ^ l^j _ where R and Ta = are respec¬ 

tively the Reynolds number and the Taylor number on scale Aq. Further 
we will consider the Reynolds number as small i? <C 1 and will construct 
on this small parameter the asymptotical development. Concerning the 
parameter D, we do not choose any range of values for the moment. Let 
us examine the following formulation of the problem. We consider the ex¬ 
ternal force as being small scale and of high frequency. This force leads to 
small scale fluctuations in velocity. After averaging, these quickly oscillat¬ 
ing fluctuations vanish. Nevertheless, due to small nonlinear interactions 
in some orders of perturbation theory, nonzero terms can occur after av¬ 
eraging. This means that they are not oscillatory, that is to say, they are 
large scale. From a formal point of view, these terms are secular, i.e., they 
create the conditions for the solvability of a large scale asymptotic devel¬ 
opment. So the purpose of this paper is to find and study the solvability 
equations, i.e., the equations for large scale perturbations. Let us denote 
the small scale variables by xo = (a;o,to) , and the large scale ones by 
X = {X,T). The small scale partial derivative operation and 

the large scale ones are written, respectively, as di, dt, Xi and dr- 

To construct a multi-scale asymptotic development we follow the method 
which is proposed in m- 
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3 The multi-scale asymptotic development 

Let us search for the solution to equations and @ in the following 
form: 


V{x,t) = ^W-i{X) + vo{xo) + Rvi + B?V 2 + H-, (4) 

H 

T(x, t) = 4r-i(X) + To{xo) + RTi + R^Ti + R^Ts + ■■■ , (5) 

H 

P{x,t) = + ^P- 2 {X) + + Po(a:o)+ 

+ R{Pi+Pi(X)) + R^P2 + R?’Pz + --- . (6) 

Let us introduce the following equalities: X = R?xo and T = R'^to which 
lead to the expression for the space and time derivatives: 


|T = a + «^v„ 

^ = a + R^dr, 




= dij + 2R^djXj + R'^dj 


dxWx^ "" tojj. 

Using indicial notation, the system of equation can be written as 

{dt + R^dT)V' + R{dj + R?Vj){V'V^) + D^EijkV'^ = 

= {dj + R^Vj)P + {dj, + 2R^diVj + R'^Vjj)V' + F^, 


(7) 

( 8 ) 
(9) 


( 10 ) 


dtT - djjT = -U" - Rdj {y*T^ , (11) 

(ft + PVi) U'= 0. (12) 

Substituting these expressions into the initial equations m and m and 
then gathering together the terms of the same order, we obtain the equa¬ 
tions of the multi-scale asymptotic development and write down the ob¬ 
tained equations up to order R^ inclusive. In the order R~^ there is only 
the equation 


ftP_3 = 0, p_3 = P-3(X). (13) 

In order R~^ we have the equation 

ftP_2=0,^P-2=P-2(X). (14) 

In order R~^ we get a system of equations: 

ftlUii - djjWU + D^eijkW\ = -(ftP-i -f ViP-3) - djWUwU, (15) 


4 



( 16 ) 


diWU = 0. 

The system of equations (O and (Hii) gives the secular terms 


- ViP_3 = 

which corresponds to a geostrophic equilibrum equation. 

In zero order we have the following system of equations: 

dtvo - djjVo + dj(WliV^o + uoTTii) + D^SijkVo = 

= -{diPo + ViP-2) + Fi, 


(17) 


diVo = 0. 

These equations give one secular equation: 

VP -2 = 0, ^ P -2 = Const. (18) 

Let us consider the equations of the first approximation R: 

dtvl - djjvl + D^EijkVi + dj{W'L-^v{ + + vWo) = 


= -VjiWUwi^) - {diPi + ViP-i), (19) 

diVi + ViWU = 0. (20) 

From this system of equations there follows the secular equations: 

ViWU = 0, (21) 

VjiWiiWi,) = -ViP-1, (22) 

The secular equations (Hzl) and dSSl) are satisfied by choosing the following 
geometry for the velocity field ( Beltrami field): 

W = (23) 

VP -1 = 0, =4> P -1 = Const. 

In the second order , we obtain the equations 


dtv\ — djjvl — 2djVjVo + dj{Wlivl + + Vov{ + vWq) 

+ D^etjkv^ = -VjiWUvi + viwi,) - {diP 2 + ViPo), (24) 

diV2 + ViVo = 0. (25) 

It is easy to see that there are no secular terms in this order. 

Let us come now to the most important order R^. In this order we 
obtain the equations 
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dtvl+dTWl^-[djjvl+2djV jv\+V j{Wl^vlWiWl^Wovl)+ 

(26) 


+dj{WU4+viwL^+vivi + v\vi + v\v{) + Eijkvl = -(ftPg + ViPi), 


diV3 + ViWi = 0 . 

From this we get the main secular equation: 

drWU - AVFli + = -ViPi, (27) 

There is also an equation to find the pressure P- 3 : 

-ViP-3 = D^eijkW!:i. (28) 

4 The velocity field in zero approxima¬ 
tion 

It is clear that the most important is equation (136 II . In order to obtain 
the se equ ations in closed form, we need to calculate the Reynolds stresses 
Vfe(uQUQ). First of all we have to calculate the fields of zero approximation 
Vq. From the asymptotic development in zero order we have 

dtvl - djjVo + W'lidkVQ + D^EijkVo = -diPo + Fq, (29) 
Let us introduce the operator Do'. 

Do = dt-djj + W'‘dk. (30) 

Using Do, we rewrite Equations (I29II : 

DqVq + D^ SijkVQ = —diPo + Pq, (^1) 

Pressure Pq can be found from condition divV = 0. 



Let us introduce designations for operators Dq: 



and for velocities: Uq = uq, Vq = vo, v§ = wq- Then excluding pressure 
from (ED) , we obtain the system of equations to find the velocity field of 
zero approximation : 

(^Do + Pxx^ uo + {^Pyx — (j^zx + ^0 = Pq 5 
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(j^xy + Dz^ Uo + (^Do + Pyy^ ^0 + (^^zy Wo = F^, 

(^Pxz — Dy^ Uo + (yPyz + Dz^ Vo + (^Do + Pzz^ Wo = Fq . ( 34 ) 

For simplicity, we choose the systeme of coordinates so that the axis Z 
coincides with the direction of angular velocity of rotation fl. Then Dx = 
0, Dy = 0, Dz = D. In order to solve this system of equations we have to 
set the force in the explicit form. Let us choose now the external force in 
the rotating system of coordinates in the following form: 

Fo = 0, Fox = /o (icos (f)2 +JCOS 01^ ; 

4>1 = klX — Wot, (f>2 = k2X — Wot, 

ki = ko (1, 0,1) , k 2 = ko (0,1,1). 

It is obvious that divergence of this force us equal to zero. Thus, external 
force is given in plane {x,y), orthogonal to rotation axis. 

The solution for equations system m can be found easily in accor¬ 
dance with Cramer’s Rule: 



Here A-is the determinant of the system m-- 


(35) 



Dq + Pxx Pyx — D 

P 

-L ZX 

A = 

Pxy D Do Pyy Pzy 



Pxz Pyz 

Do -\- Pz 



FS Pyx -p 

Pzx 

Ai 

= 

Fq Po + Pyy 




0 pyz 

Do + Pzz 



Do -|- Pxx Fq 

p 

j- ZX 

A2 


Pxy + P Fq 




Pxz 0 

Dq + Pzz 



Po T Pxx Pyx 

-D FS 

A3 

= 

Pxy H" D Dq + Pyy Fq 



Pxz Pyz 0 


(36) 


(37) 


(38) 


(39) 


After writing down the determinants in the explicit form, we obtain: 


no = i [( 5 o + Pyy) [Do + Pzz) - {Pyz) (p.,)] FS + 

+ i [(p..) (p..) - (p.. - d) {Do + P..)] Po^ 

no = i [(P.3) (P.^) - (P., + P) (Po + P..)] Po"+ 
+ i [(Po +P..) (Po +Pzz) - {Pxz) (Pz.)] Po^ 
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w;o = ^ + d) (Py.) - (p,.) ( 5 o + Pyy)] FS+ 

+ s [(p-) {Py- - Pj - (Po + P-j (p«-)] Pq . 

A = ^Do + Pxx\ I^^Po P Pyy^ ^Po P Pzz^ — (^Pyz^ — 

— (Pyx — DJ j^^Pr;/ P P^ ^Po P Pzz^ “ (^Pxz'j P (^3) 

P (p^a^) [(Pa:B P P) (Pj'^) “ (Po P Pj/a) (P^^)] ■ 

In order to calculate the expressions (|40ll-([43ll we present the external 
force in complex form: 


Fo" = ^ , Pd' = Y ■ (^4) 

Then all operators in formulae (I40I) - (I42I) act from the left on their eigen¬ 
functions. In particular: 

i5oeP= = eP=Po (k2, -ojo) , Poe^i = eP^Po ffci, -wo) , 

\ ^ X ^ (45) 

Ae*<#>2 = eP2 A , -cjo j , A (^ki , -wo j . 

To simplify the formulae, let us choose fco = 1, wo = 1. 

We will designate 

Do {k2, —wo^ = 2 + i {wy — 1) = Ay, 

Do {ki, -woj = 2 -I- i {wx - 1) = Ax. (46) 

Before do further calculations, we h ave to note that some components of 
tensors Pij {k^^ and Pij {k 2 ^ vanish. Let us write the non-zero compo¬ 
nents only: 

Pyx (ki^ = ^P, Pxz (^k2^ = --j^D, Pxy {k'^ = -^P’ 

P.. (fci) = \d. (47) 

Taking into account the formulae (|4^-(|47ll. we can find the determinant: 


(t.) 


— Ax P 2P 


( 6 ) 


In a similar way we find velocity field of zero approximation: 

p«02 £) g«01 


r Ay 

Wo — fo 


2 Al + lD^ 


Vo = - fo 


D 


p»<A2 


4 Al + 1P2 


P/o 


P/o- 


4 42 P ip2 




2 Ai + ip2 


-f C.C., 
-f C.C., 


UIQ = fo 


D 


o*</>2 


4 Al + ip2 


-/o 


P 




4 Ai + ip2 


-f C.C.. 


(48) 


(49) 

(50) 

(51) 
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5 Reynolds stress and large scale insta¬ 
bility 


To close the equations 123 we have to calculate the Reynolds stresses 
wqUo and wovo ■ 

These terms are easily calculated with help of formulae (gSll-ISI]). As 
a result we obtain: 


il 

2 \Al + lD^ 


/o 

8 |Ai + iZ32 


/o 


8 \Al + 


fo D 
2 \Al + \D^ 


Now equations (|27l) are closed and take form: 


OtWj: — AWj; + —WqUo = 0, 
oz 


(52) 


drWy — AWy — -^wovo = 0. (53) 

oz 

We calculate the modules and write the Reynolds stresses ()52l) in the 
explicit form: 


WqUq = 


fl _ 

2 16 (Wy 


_ D _ 

lf+[4+^D^-{Wy 


1)1 


2 


il __ 

8 I6(w,-1)2+[4 + iD2-(w,-l)2]2’ 

__fl _ il _ 

^ lQ{wy-lf+[A+\D^-{wy-lfY 


il _ D _ 

2 I6(w, - 1)2 + [4+- (w, - l)2]2' 


(54) 


With small W^, Wy Reynolds stresses dSH) can be expanded in a series in 
the small parametersWc, Wy. Taking into account the formula: 


\Wi,y + hD'^\ 


= Const — 


32 (1)2 - 10) 
[(7)2+6)2 + 64]' 


we obtain the linearized equations (1531) : 


d 


^2 


-w-_ 

dT dz^ 

d 92 


q/qD d 
2 dz 

q/o7)2 9 
8 dz 


Wy + 


Wy + 


afoD^ d 

8 dz 

afoD 9 
2 dz 


ir. = 0 , 


W^ = 0. 


32 (10 - 7)2) 

[(7)2+6)2 + 64]^' 

We will search for the solution of linear system (1551) in the form: 


(55) 
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Figure 1: Phase portrait of the dynamical system ( 1551 ) . with D = 1, Ci = —0.03, 
C 2 = 0.03. One can see two hyperbolic singular points and stable and instable 
knots. 


Wx,Wy ~ exp (yT + ikZ). (56) 

We substitute (I56II in equation (I55II and obtain the dispersion equation: 


, = -ik^±k^-k\ 


(57) 


The dispersion equation dsa shows that equation system ([55]) has instable 
oscillatory solutions with oscillatory frequency ui = fe— g— and instability 
growth rate 7 = — k^. The instability is large scale because the 

instable term dominates over dissipation on large scales: — > k. The 

0,2 ^ 4^2 

maximum growth rate of instability is equal to 7 max = —^—, and is 

achieved on the wave vector fcmax = —. As a result of the development 

of instability the large scale helical circular polarized vortices of Beltrami 
type are generated in the system. 


6 Saturation of instability and nonlinear 
vortex structures 

It is clear that with increasing of amplitude nonlinear terms decrease and 
instability becomes saturated. Consequently stationary nonlinear vortex 
structures are formed. To find these structures let us choose for equations 
dSlI ^ = 0 and integrate equations one time over Z. We obtain the 
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Figure 2: The kink which connects the hyperbolic point with stable knot with 
= 1, Cl = 0.04, C 2 = 0.04. When approaching the stable knot one can see 
rotations of velocity field. 


system of equations: 


— WT = wqUo + Cl, 
^ Wy = wm+C2. 


From equations (I58II follows: 


dwx _ wouo + Cl 
dwy Wo Vo + C 2 ’ 

After integrating the system of equations in we obtain: 


(58) 


(59) 


j wovodwx + C 2 WX = j wouodwy + CiWy. (60) 

Integrals in expression (IMI are calculated in elementary functions (see 
[H]), which give the expression for first integral of motion J of equations 

ISSl): 

J _ _ 

8 [4+ 

D (wx — 1)^ + (wx — 1) D 1/2 + 4 + 

+ 25/2 (8 + ^2) In (^^_l)2_(^^_i)^y2 + 4+i02 + 

D 1)2 _ 4.1^52 

- 
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Figure 3: The kink which connects the instable and stable focuses with D = 1, 
Cl = 0.04, C 2 = 0.04. One can see the internal helical structure of the kink. 




+ 


D (Wy - 1) + {Wy — 1) D\/2 + 4 + 

■ In -^^-;=-T—- + 


25/2 (8 JJ 2 ^ ^^2 + 4 + iD 2 


+ 7 


D 


■ arctan 


(^^_ 1 ) 2 _ 4 _ 1^2 


+ ClWy + C2WX 


8(8 + L>2) 4(wy-l) 

Equations (I58II can be easily calculated numerically using standard tools. 
In particular, this allows to construct phase portrait of the dynamical 
system (I58II ( figH]) and to get the most interesting solutions which link 
singular points on phase plane. See for example figlU where the hyper¬ 
bolic singular point is connected with the stable knot and fig.3, where 
the solution connects instable and stable focuses. All these solutions cor¬ 
respond to the large scale localized vortex structures of kink type with 
rotation, generated by the instability which has been found in this work. 


7 Conclusions and discussion of the re¬ 
sults 

In this work we find new large scale instability in rotating fluid. It is 
supposed that the small scale vortex external force in rotating coordinates 
system acts on fluid which maintains the small velocity field fluctuations 
(small scale turbulence with small Reynolds number R, R 1 ). For the 
real applications this Reynolds number should be calculated with help of 
the turbulent viscosity. The asymptotic development of motion equations 
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by small Reynolds number allows obtaining motion equations for the large 
scale. These equations are of the hydrodynamic a-effect type, in which 
velocity components WT, Wy are connected by the positive feedback. This 
may result in the appearance of the large scale vortex instability. The large 
scale vortices of Beltrami type are formed due to this instability in rotating 
fluid with small scale exterior force. With further increase of amplitude 
the instability stabilizes and passes to a stationary mode. In this mode 
the nonlinear stationary vortex structures form. Different vortex kinks 
belong to the most interesting structures. These kinks link stationary 
points of dynamical system (1581) . The kink which links hyperbolic point 
with stable knot has rotations around the stable knot as shown on figl21 In 
the kink which links instable and stable focuses, vector field turns around 
two singular points, see figO 

Let us note that unlike previous works about hydrodynamic a- effect in 
rotating fluid, the use of the asymptotic development allows to construct 
naturally the nonlinear theory and to study the stationary nonlinear vor¬ 
tex kinks. 
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